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ABSTRACT 



Context. The detection and identification of oscillation modes (in terms of their f, m and successive n) is a great challenge for present and 
future asteroseismic space missions. The "peak tagging" is an important step in the analysis of these data to provide estimations of stellar 
oscillation mode parameters, i.e., frequencies, rotation rates, and further studies on the stellar structure. 

^ Aims. To increase the signal-to-noise ratio of the asteroseismic spectra computed from time series representative of MOST and CoRoT 

O observations (30- and 150-day observations). 

Methods. We apply the curvelet transform - a recent image processing technique which looks for curved patterns - to echelle diagrams built 
using asteroseismic power spectra. In this diagram the eigenfrequencies appear as smooth continuous ridges. To test the method we use Monte 
Carlo simulations of several sun-like stars with different combinations of rotation rates, rotation-axis inclination and signal-to-noise ratios. 
Results. The filtered diagrams enhance the contrast between the ridges of the modes and the background allowing a better tagging of the 
modes and a better extraction of some stellar parameters. Monte Carlo simulations have also shown that the region where modes can be 
detected is enlarged at lower and higher frequencies compared to the raw spectra. Even more, the extraction of the mean rotational splitting 
from modes at low frequency can be done more easily than using the raw spectrum. 
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1. Introduction tometry missions and ground-based velocity networks must 

be used to provide observations of stellar oscillations with- 

HeUoseismology - the study of solar oscillations - is a ^^^ ^^^^^ limitations. With the current MOST' and WIRE^ 

powerful probe of the structure and dynamics of the Sun ^^^^jj^^^^ ^^^ ^^^ ^^^^^.^ (.qrq^3 ^^^^^^^ asteroseismology is 

which has provided great improve ments in our understand - blooming. However, we still have to deal - in the case of solar- 

ins of stellar evolution and st ructure < Turck-ChiezeetaHl993t ^^^ oscillations - with very small signal-to-noise ratio (here- 

I Christensen-D alsgaarcf 2002, and references therein). Those ^^^^^. ^/^^ observations as a consequence of the weakness of 

successes push the community to apply seismic techniques to ., , • •. ... . , . . , ,■ ,, 

^ .; Fi' J M j-jjg luminosity variations. Moreover, stars cannot be spatially 

other stars, opening the doors to asteroseismology, the study ^^^^j^^^ ^^^ ^^^^ gj^^^j oscillation modes can be observed, 

of stellar oscillations. These oscillations have akeady been ob- ^^ ^^^.^.^^^ ^^ ^^^^^^ ^^^^ ^^^^^^ ^^ ^^^ ^^^^^.^^ ^^^^^ ^^^ 

served from ground and space. The ground-based observations ^^^ notation-axis inclination separately. Without knowing these 

are limited by the day-night cycle, which introduces aliases ^^^ j^^^ ^^^jj^^ properties, the tagging of the modes in terms 

in the observations, but allow to use Doppler velocity mea- ^ -.u • »• /^ \ j • u » i 

^^ ^ of their properties (t, m) and successive n may be extremely 

surements. They have provided dat a with sufficient quality ^jg^^^j^ ^^ ^^^^^ ^^^ ^^^ ^^.^^^^^ ^^ ^^^^ ^^jj ^^^ ^^ ^^ ^^ 

to detect solar-like oscilla tions (see | Bouchy & Carriei| | 2003t ^j^^ p^^,^^ ("peak-bagging") but to provide a good description 
iBedding & Kieldserill20 03' and references therein). To reduce 

the aliases, multi-site campaigns have been carried out but they 

are too short to have a good frequency resolution. Space pho- ' Microvariability and Oscillatio ns of STa rs i'MatthewJl99a) 

2 Wide-field Infra Red Explorer JSuzasTet aU2000,) 

Send ojfprint requests to: P. Lambert ' Convection Rotation and planetary Transits jBaglin et all200ll) 
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of the model to be fitted after having put the correct labels on 
the modes ("peak tagging"). To do this, it has been proposed 
to use the echelle diagram where the modes follow ridges de- 
pending on the stell ar properties. To improve the S/N ratio 
bedding et al.l ( 120041) proposed to filter this diagram by a verti- 
cal smoothing. However the smoothing works well only when 
the ridges are quasi-vertical which means a very good a pri- 
ori knowledge of the large difference and is restricted to the 
asymptotic part of the spectrum. We propose here to follow a 
similar approach but using new mathematical denoising tech- 
niques better suited to the study of curved ridges. 

At the end of the last decade, the application of mathemat- 
ical transforms based on wavelets to analyze astronomical im- 
ages has been widely developed. The first wavelet algorithms 
were well adapted to treat images with isotropic elements. 
However, this description presented a limitation in the context 
of astrophysics, where objects such as filaments or spiral struc- 
tures exhibit a highly anisotropic character (in shape and scale). 
New tr ansforms, the ridgelet tCandes_1998.) and cur velet trans- 
forms dCandes & Donoholl999HStarck et alJ2002l) . were then 
developed to deal efficiently with such objects. Astrophysical 
applicati ons (image denoising) of t his technique have been pre- 
sented in lStarck et alJ (I2003ll2004l) to analyze images of gravi- 
tational arcs, the Saturn rings or the CMB (Cosmic Microwave 
Background) map. 

In this paper we suggest to use the curvelet transform to 
analyze asteroseismic observations (more precisely the stellar 
echelle diagrams), in order to improve the "peak tagging" of 
the oscillation modes and even the resultant "peak bagging". To 
illustrate the application of this denoising technique in the as- 
teroseismic case, we have performed Monte Carlo simulations 
of ideal asteroseismic data contaminated by different levels of 
stochastic noise. We start in Sect. 2 by a quick reminder of the 
properties of stellar oscillation modes in the solar-like case and 
the construction of the echelle diagram. In Sect. 3 we introduce 
multiscale transforms, in particular the ridgelet and the curvelet 
transforms. In Sect. 4, the simulated data of a star with an os- 
cillation spectrum similar to the Sun but with different rotation 
axis inclinations and rotation rates, are presented. In Sect. 5 we 
discuss the results obtained in the simulations. 



2. Properties of solar-like oscillations 

Only low-degree stellar oscillation modes can be detected 
and observed with the present generation of instruments. The 
asymptotic theory of oscillation modes (n » €) is then ad- 
equat e and can be used to study them. Fi rst order (iTassoull 
1980l) and second order developments IV^rontsovl Il99l| 



Lopes & Turck-Ch ieze 1994; Roxburgh & Vorontsov 2000a'b' 



have been made to describe solar and stellar oscillations. In the 
case of solar-like stars, where p-modes predominate, the fre- 
quencies can be developed as: 
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Fig. 1. Portion of the theoretical spectrum (top) and echelle di- 
agram (bottom) for a sun spinning ten times faster than the Sun 
and seen under an angle of 50°. This is the ideal power spec- 
trum used in the simulations described in Sect. 5. 

in this expression t and n are respectively the degree and the 
radial order of the modes and 
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c, is the internal stellar sound speed, a is a phase-shift term 
and i/' is a function which allows to tak e into account the grav- 
itation al potential in the central region (ILopes & Turck-Chiezel 
I1994I) . From the asymptotic approach, we can extract general 
properties of modes and better understand the physics hidden 
in the frequencies behavior. The large frequency spacing, de- 
fined as Av„/ - v„+i/ - v„x, tends asymptotically to Avq, re- 
lated to the mass and radius of the star; the small frequency 
spacing, (5f/+2V - v,,/ - Vn-iz+i, can be approximated to first 
order by i4( + 6)Avo/(47:\„j) J^* ^f . This variable is re- 
lated to the derivative of the sound speed and enhances the ef- 
fect coming from the central regions, providing constraints on 
the age of the star. Finally the second diff'erence is defined as 
^2V = Vh+i/ - 2v„/ + Vn^\,t- Its Variation s provide information 
about the extent of the convective zone JMonteiro et al.ll2000t 
Ballot et al. 2004b) or the helium abundance in the stellar en- 
velope CBasu et al_2004) . 

Under the rotation effects the azimuthal order m (-£ < m < 
{) is needed to characterize the oscillation spectrum. If the an- 
gular velocity Q. is uniform (Ledoux 1951) . the mode frequen- 
cies are asymptotically approximated by: 
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V«,f,K 



v„/ + mD./27T — v„_f + mdv 



(2) 



p. Lambert et al.: Curvelet analysis of asteroseismic data 



where 6v is the rotational splitting. Equation |2 shows that 
modes are (2( + l)-times degenerated among the azimuthal 
order: a single peak in the spectrum becomes a multiplet. Its 
corresponding structure depends on the rotation rate, the incli- 
nation axis of the star and its stochastic excitation. The solar- 
lUce mode lifetimes (a few days) are expected to be much 
shorter than the length of the future space observations (a few 
months). In consequence, the relative amplitude ratios inside 
a multiplet will only depend, in average, on the inclination 
angle and the spacing between these different m-components 
( iGizon & Solankill2003h . Thus if the different m-components 
of a multiplet can be identified and tagged with the correct 
{{, m), they can provide a good estimation of both the rotation- 
axis inclination / and the rotational splitting 6v, allowing a bet- 
ter mode parameter extraction through the fitting of the spec- 
tra. The effect of the stochastic excitation on an isolated mode 
could be minimized by computing the average of these param- 
eters on several modes (see for example the n-collapsogramme; 
[Ballot et al. 2004a). 

Equation n shows that the even {( - 0,2) and odd {{ - 
1 , 3) modes have respectively almost the same frequency, only 
separated by the small spacing di^+iv. In addition, they are 
separated regularly in frequency by the large spacing Av„/. 
This property allows us to build the so-called echelle diagram 
( Grec et al. 1983), which is currently used to identify modes for 
solar-like oscillations. It is a 2D representation of the power 
spectrum where this one is folded onto itself in units of the 
large spacing. In such representation the modes appear as al- 
most locally vertical ridges (see Fig. ^. The echelle diagram 
is a powerful tool for the "peak tagging" since assigning the 
correct {t, m) values to the peaks is easier when the multiplet 
structure is well identified in this diagram. The successive n 
values are obtained from each individual horizontal line. 

3. Multiscale Transforms 

3.1. The Wavelet Transform 

The wavelet transform provides a framework for decomposing 
images into their elementary constituents across scales by re- 
ducing the number of significant coefficients necessary to rep- 
resent an image. The continuous wavelet transform of a 2D sig- 
nal is defined as: 



W{a,b„bj) = -^ JJ/(x,3')r {^^ 



■hi y-bj 



dxdy 



(3) 



where W{a, b) are the wavelet coefficients of the function fix), 
ij/(xy is the conjugate of the analyzing wavelet, a > is the 
scale parameter and b is the position parameter The continu- 
ous wavelet transform is the sum over all the positions of the 
signal fix, y) multiplied by the scaled and shifted versions of 
the wavelet i/'((x - bi)/a,iy - bj)/a) (cf. Fig. |2] top panels). 
This process produces wavelet coefficients that are a function 
of scale and position. 

However, the classical wavelet transform only address a 
portion of the whole range of interesting phenomena: isotropic 
features at all scales and locations. One of the drawbacks of the 






Fig. 2. Examples of 2D wavelets (top panels) and ridgelets (bot- 
tom panels). The top right wavelet has a greater scale parameter 
than this on the left. The bottom right ridgelet has different ori- 
entation and width than the left one. 



two-dimensional wavelet transform is that it does not achieve 
an efficient analysis of images which present high anisotropy. 
For instance, the wavelet transform does not efficiently approx- 
imate 2D edges, since a large number of large wavelet coeffi- 
cients, scale after scale, are required, making difficult its analy- 
sis. In order to solve this problem two new mathematical trans- 
forms, namely the ridgelet transform and the curvelet trans- 
form, were introduced. 



3.2. The Ridgelet transform 

The ridgelet transform was developed to process images in- 
cluding ridges elements ICandesll998l) . It provides a represen- 
tation of perfectly straight edges. Given a function fixi,X2), 
the representation of this latter is the superposition of elements 
of the form a^^^^i/zHxi cos 6 + X2 sin0 - b)/a), where (A is a 
wavelet, a > a scale parameter, b a location parameter and 
an orientation parameter. The ridgelet is constant along lines 
xi cos 6 + X2sin6 - const, and transverse to these ridges it is 
a wavelet. Thus, contrary to a unique wavelet transform, the 
ridgelet has two supplementary characteristics: a length, equal 
to this of the image and an orientation, allowing the analysis of 
an image in every direction and so exhibiting the edge struc- 
ture. Fig. 121 (bottom panels) shows two examples of ridgelets. 
The problem is that in the nature edges are typically curved 
rather than straight so ridgelets alone cannot yield an efficient 
representation. 



3.3. The Curvelet transform 

3.3.1. Description 

Ridgelets can be adapted to represent objects with curved edges 
using an appropriate multiscale localization: at a sufficiently 
fine scale a curved edges can be considered as almost straight. 
Candes & Donoho ( 1999) developed the curvelet transform us- 
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Image 




Fig. 3. Sketch illustrating the curvelet transform applied to an 
image. The image is decomposed into subbands followed by 
a spatial partitioning of each subband. The ridgelet transform 
is applied to each block. The finest details correspond to the 
highest frequencies. 



ing ridgelets in this localized manner. Fig.|3shows the different 
steps of the curvelet analysis of an image: 

1. Image decomposition into subbands: as a set of wavelets 
bands through a 2D isotropic wavelet transform. Each band 
corresponds to a different scale. 

2. Smooth partitioning: each subband is partitioned into 
squares - blocks -, whose size is appropriate to each scale. 
The finest is the scale, the smaller are the blocks. 

3. Ridgelet analysis: it's applied to each square. 

The implementation of the curvelet transform offers an ex- 
act reconstruction and a low computational complexity. Like 
ridgelets, curvelets occur at all scales, locations and orien- 
tations. Moreover contrary to ridgelets, which have a given 
length (the image size) and a variable width, the curvelets have 
also a variable length (the block size) and consequently a vari- 
able anisotropy. The finest the scale is, the more sensitive to the 
curvature the analysis is. As a consequence, curved singulari- 
ties can be well approximated with very few coefficients. 

3.3.2. Denoising images: filtering curvelet coefficients 

To remove noise a simple thresholding of the curvelet coeffi- 
cients has been applied to select only significant coefficients. 
One possible thresholding of a noisy image consists in setting 
to all non-significant curvelet coefficients cijj, i, j and I re- 
spectively the indexes of the line, row and scale: it is the so- 
called hard-thresholding: 



CiJJ 



1 if Cijj is significant 
if Cijj is not significant 



(4) 



Commonly, c, ^ / is significant if the probability that the curvelet 
coefficient is due to noise is small, i.e., if the curvelet coefficient 
is greater than a given threshold. A basic problem remains: the 
choice of the threshold. Usually, this threshold is taken equal 



to kcTj, where ctj is the noise standard deviation at the scale j 
and fc is a constant taken equal to 5 in our filterings. 

Simple thresholding of the curvelet coefficients is very 
competitive ( Starck et al. 2002) with "state of the art" tech- 
niques based on wavelets, including thresholding of decimated 
or undecimated wavelet transforms. 



4. Simulation of data 

To characterize the curvelet denoising technique applied to the 
asteroseismic data, we have simulated typical solar-like obser- 
vations varying different parameters: S/N ratios, observational 
lengths, rotation-axis inclinations, rotation rates... With this ap- 
proach we know the input parameters in advance and we can 
evaluate the quality of the results given by the curvelet analysis 
and its limits. 

In the simulations shown in this paper, we use the oscilla- 
tion spectrum of a star similar to the Sun but seen under differ- 
ent conditions. Different rotation-axis inclinations (/ = 50° and 
90°) and rotation rates {Q. = Q.q, SQq, and IOO0) have been 
considered. An ideal power spectrum were constructed first. 
Only the modes f < 3, n = 12-25 were simulated. The mode 
parameters - frequencies (v), amplitudes (A) and widths (F) - 
were obtained from the a nalysis of GOLF (G lobal Oscillations 
at Low Frequency) data (Garcia et al. 2004). The amplitudes 
were corrected to take into account the difference between in- 
tensity and velocity observations. Modes were simulated with 
symmetrical Lorentzian profiles as the asymmetry is expected 
to be at the level of the noise. Following the method described 
in lFierrv Fraillon et alJlll998l) . a multiplicative noise, sl^^ with 
2 d.o.f. statistics, has been introduced to reproduce the stochas- 
tic excitation of such modes (see also Anderson et al. 1990) . 
The S/N ratio of the "resultant" raw power spectrum was de- 
fined as the maximum of the bell-shaped p-mode power (i.e. 
the highest simulated p mode) divided by the noise dispersion. 
The simulated background is flat assuming that it has been pre- 
viously fitted and removed as it is usually done for the Sun 
faarvev 1985) . 

Several Monte Carlo simulations have been performed for 
each ideal spectrum. Realistic S/N, with values ranging from 
5 to 15, have been used to cover a wid e range of situation s 
(compatible with what it is expected, (see Bariin et alJ200ll) '). 
In each realization of the Monte Carlo simulation the same 
level of noise has been randomly added to the correspond- 
ing ideal spectra. Therefore all the realizations, in a given 
Monte Carlo simulation, have the same S/N ratio. The simu- 
lated spectra have been computed for two resolutions, « 0.38 
and ^ 0.077 juHz, corresponding respectively to 30-day and 
150-day observations. The first are representative of MOST ob- 
servations and the short CoRoT runs while the latter are of the 
same length than the long CoRoT runs. 

Simulations of other stars, like some potential main CoRoT 
targets, with different masses, ages and, in consequence, inter- 
nal structures have been made. The results have already been 
presented and discussed during the CoROT workshops #8 and 
#9 obtaining the same qualitative results. For the sake of clarity, 
they are not shown here. 
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Fig. 4. Effect of the curvelet denoising on the mode visibility for S /N - 5. Each picture shows 120 realizations out of the 500 
done in our Monte Carlo simulation. Each horizontal line corresponds to a single realization. The top panel is the raw spectra 
and the bottom is the curvelet filtered one. 



5. Discussion 

Once the spectra have been computed, the echelle diagrams can 
be built with a fixed folding frequency. This one corresponds to 
the mean large frequency spacing Avq, identified either by com- 
puting the EFT, the autoco rrelation of the spectra or any other 
technique (see for example Regulo & Roca Co rtes 2002). The 
denoising based on the curvelet transform is then applied to this 
echelle diagrams. It is important to note that artifacts may ap- 
pear in the filtered spectra at frequencies v*=vo-i-A:Av(), with k an 
integer, when random small structures appear in the echelle dia- 
grams. However, their appearance and position strongly depend 
on the folding frequency and are very sensitive to its value. 
Therefore they can be easily identified. The artifacts can be re- 
duced (in contrast to the regions containing signal) by building 
echelle diagrams with slightly different folding frequencies and 
averaging the resultant filtered spectra. 

In order to present the results of data analysis using the 
curvelet denoising method, we have selected the case of a sun- 
like star seen with an inclination angle / - 50° and with a rota- 
tion Q = 1OQ0. A portion of the ideal spectra constructed for 
this star can be seen in Fig.n(top panel). Monte Carlo simu- 
lations were then performed, giving rise to different sets (each 
one with 500 realizations) of raw spectra with different S /N 
ratios. The echelle diagrams were constructed using a folding 
frequency of 135.18 yuHz, obtained by computing the EFT of 
the raw spectrum. 



5.1. Peak tagging 

In those cases, with a high S /N (typically 15), the mode struc- 
ture is clearly visible in each raw spectrum and also on the 
echelle diagram. The different ridges can be easily identified 
and tagged. Although the filtering gives enhanced denoised dia- 
grams and unfolded spectra, it does not contribute significantly 
to the mode identification. 

In the lower S /N cases, however, the situation is different. 
Figure |4] shows some of the results of the Monte Carlo sim- 
ulation for S /N=5. The upper panel corresponds to 120 real- 
izations among the 500 computed for the raw spectra in the 
frequency range 2450-2920 //Hz. Each horizontal line corre- 
sponds to a single realization. Some patterns can hardly be 
seen. The lower panel represents the same spectra after apply- 
ing the curvelet filtering. A series of vertical ridges clearly ap- 
pears. From the left to the right on the panels, they can be iden- 
tified as the ({ = 2; m = ±1), the ^ = (blended with the ^ = 2; 
m - +2) and the {{ - I; m = -1,0,-1-1). The improvement 
of the contrast is important in all the realizations and allows to 
distinguish the different components of a mode, making easier 
the identification and the tagging. 

The identification is harder when looking at each individ- 
ual spectrum and requires the use of the echelle diagram. Fig.|5] 
shows an example of raw (left) and filtered (right) 150-day ob- 
servation power spectra (top and middle panels) and the cor- 
responding echelle diagrams (bottom panels) for a S /N - 5 
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Fig. 5. Raw (left) and filtered (right) power spectra (top and middle panels) and echelle diagrams (bottom panels) for aS /N - 5 
realization. The short dashed lines in the power spectra represent the position of the theoretical frequencies. From left to right, the 
three first equidistant lines indicate the components ;« = -1, 0, 1 of /" = 1 modes, the two next indicate the strongest components 
of ( = 2 (m - -1 and 1), and the the last indicates ^ = 0. In this case only two components of the { = I and the /" = mode are 
slightly visible in the raw diagram. On the curvelet filtered one, the three ^ = 1 components appear as well as the { = Q and the 
components m - ±1 of the { - 2 modes. 



realization. Input frequencies are indicated by the short dashed 
lines above the spectra. The mode peaks can hardly be distin- 
guished in the raw spectrum and can easily be confused with 
noise. For the range 2780-2920 //Hz, only a strong peak at 
2900 /iHz can be considered not to be noise. In the region 
3060-3180 fjHz the peaks are visible and we can attempt to 
identify the £= 1 and {-0 modes but it is still unclear On the 
contrary, on the corresponding parts of the filtered spectrum, 
the structures of the {= 1 mode with three components, the {=0 
mode and even the strongest components of the {=2 mode are 
visible. The raw echelle diagram gives no extra information be- 
cause of the very weak ridges and low contrast with the back- 
ground. The weakest components can hardly be detected and 
no tagging can be done. The curvelet filtering provides a con- 
trast enhancement of the ridges on the echelle diagram. Thus 
three almost equidistant strong ridges appear on the left of 
the diagram and one strong ridge with two weaker ones on 
the right. The corresponding patterns can be seen on the fil- 
tered spectrum corresponding well to the theoretical frequen- 
cies. Since the modes { = 3 are not visible, and according to 
the amplitude of the strongest peak on the left, we can suggest 
that the three strongest peaks correspond to a ^ = 1 multiplet 
and the other ones to the /' = 2 and { - Q modes. 



Consequently, when the tagging is done it is also easier to 
have a first estimation of both the mean rotational splitting and 
the rotation-axis inclination, since the visibility of the multi- 
plet is increased. From the spacing between the components of 
the mode C - \,a. first estimation of the mean rotational split- 
ting of the star can be done, as well as an estimation of the 
inclination angle, according to their relative amplitude ratios. 
We have selected the extraction of one parameter: the mean 
rotational splitting of the i- 1 mode at low frequency (2540- 
2550 juHz), to quantify the improvement of the curvelet fil- 
tering. This region is particularly interesting because the line 
width is still small and the modes, when they are visible, can 
be easily identified. Thus, in a sample of 100 realizations of 
the Monte Carlo simulation, we have obtained in 90 of them 
a better estimation of this parameter in the filtered spectra. In 
fact, in the raw spectra it was very exceptional to obtain a good 
result. With the filtered spectra a mean rotational splitting of 
{6v) - 4.05 ± 0.30 //Hz was found, which is very close to the 
actual splitting included in the ideal spectra {6v) - 4.0 //Hz. In 
addition, specific methods can be applied to improve the extrac- 
tion of these parameters by using different strategies of spec- 
tra fitting as the on es developed by Gizon & Solanki ( 2003j ) or 
Ballot et alJ ( l2006l) . In the case of the 30-day observations, the 
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curvelet filtered echelle diagram is still very noisy and it does 
not help in recognizing the ridges. However the correspond- 
ing denoised power spectrum is much better despite the lower 
resolution (5 times less than in the long runs), even for small 
S /N ratios (~ 5). The modes /" = 0, 2 and i - 1 can be distin- 
guished, at the maximum power, while it is not obvious to do 
so in the raw spectra. Therefore, we consider that a 30-day run 
is the minimum length needed to have reliable results with the 
curvelet denoising technique. 

iGarcia et al._c2006.) analyzed the first available MOST pub- 
lic Procyon A data (32-day obs ervation) using t he cur velet 
technique. Previous analysis by Matthews et al. (2004) did 
not reveal the presence of any p-mode structure in this star 
Therefore, due to its tiny S/N ratio the results of the curvelet 
denoising should be taken with care. Nevertheless, an excess 
of power seems to appear in the region where it is expected and 
taking the 15 most prominent peaks in this region, many are in 
agreement, inside the error bars, with previous tagged modes 
using ground-based velocity observations. 



5.2. Extraction of p-mode parameters 

Once the mode identification and tagging are done, the extrac- 
tion of the mode parameters can be performed. To illustrate 
how this extraction can be improved by using the denoised 
spectrum we have extracted the central frequency of the modes 
in both the raw and the filtered spectra. To determine this pa- 
rameter, modes have been fitted by Lorentzian profiles using 
a maximum-likelihood estimator in the classical way: adjacent 
pairs of even (^ = and £ -2) modes are fitted together, while 
^ = 1 is fitted alone, due to the small amplitudes of ^ = 3 
modes. For each multiplet, the fitted parameters are the cen- 
tral frequency v„/, the amplitude A„/, the linewidth f,,/ and 
the background b. The amplitude ratios inside the multiplets 
and the rotational splittings have been fixed thanks to the pre- 
liminary estimation done in the previous section (cf. 5.1). The 
fitting procedure provides for each adjusted parameter X an as- 
sociated error cr(X) computed by Hessian-matrix inversion. 

The raw spectra follow a^^ with 2 d.o.f. statistics, whereas 
the filtered spectra have a x^ with a higher d.o.f. statistics 
(close to a Gaussian distribution depending on the number of 
filtered coefficients). According to Appourchaux (2003), it is 
possible to fit spectra following a x^ with more than 2 d.o.f. 
statistics with a classical procedure developed for ax^ with 2 
d.o.f. statistics: parameters are correctly fitted, but computed 
errors have to be adapted a posteriori. However in our case, 
due to filtering, points of filtred spectra are correlated (we have 
estimated that one point is correlated with ~10 neighbouring 
points). This correlation should have to be considered, but we 
have neglected its effect on the fitting procedure in the present 
study. This assumption is validated by the Monte Carlo simula- 
tions. Such a global filtering induces also correlations between 
the different lines of the echelle diagram. Thus the errors on 
parameters of different modes (typically (n,{) and (« + 1,^)) 
can be correlated. These correlations will have to be taken into 
account especially during the comparison of frequencies ex- 
tracted by this way to stellar models. 



From the 500 reaUzations of the Monte Carlo simulation, 
we derived for each mode and for both the raw and the fil- 
tered spectra the mean value of the extracted frequencies (v„/), 
their mean computed errors (cr(v„/)) and the dispersion of fre- 
quency distribution o-*{v„/) (the real error). We have verified 
that o-*(v„/) « (cr(v„/)) for fits performed on the raw spectra 
and we have cr*(v'„ f) < (cr(v„/)) for fits performed on the fil- 
tered ones. As expected, the error bars on the fitted frequencies, 
computed by Hessian-matrix inversion, are overestimated. 

Figure |6l shows the difference between the mean fitted fre- 
quencies (v,,/) and the theoretical frequencies v,„ of the simu- 
lated star discussed in the previous section (S/N = 5). The error 
bars correspond to the dispersion o-*{v„(). For each £, the error 
bars of the filtered spectra are smaller than those of the raw 
spectra. In addition, the range where modes can be detected, 
tagged and fitted is extended. While the difference {v„/) - v,„ 
is only flat in the central region of the raw power spectrum 
(e.g. for ^ = 0, in the range n = 18-22), it extends at higher 
and lower frequencies (e.g. for { - 0, the range is extended to 
n - 16-23) in the filtered one. 

6. Conclusions 

The application of a noise reduction technique based on the 
curvelet transform to echelle diagrams improves the identifica- 
tion - "peak tagging" - of stellar acoustic modes. In observa- 
tions with aS IN ratio as small as 5 we are still able to recover 
the mode pattern and extract reliable asteroseismic information 
in both small and long runs (30-day and 150-day observations 
respectively). Below this S/N and with shorter observations, the 
method efficiency is reduced drastically. The rotational split- 
tings and the rotation-axis inclination can be better estimated 
using the filtered spectrum. In particular, Monte Carlo simu- 
lations showed that a better extraction of the mean rotational 
splitting from modes at low frequency can be done in 90 out of 
100 realizations using the filtered spectra. The uncertainty on 
the extracted rotational splitting of a typical sun-like star seen 
with an inclination angle / - 50° and with a rotation Q = \OSIq 
is very small, ~0.30 /iHz. These parameters can then be used 
to have a set of guesses or a priori values to perform individual 
fits of the spectra. We have also shown that the range of the 
frequency extraction can be extended at higher and lower fre- 
quencies using the filtered spectra. Finally, simulations of the 
short run observations have demonstrated that this method can 
also be applied to lower resolution spectra with good results. 
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